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Abstract 

^ ' The definition of entropy obtained for stationary black holes is ex- 

tended in this paper to the case of non-stationary black holes. Entropy is 
defined as a macroscopical thermodynamical quantity which satisfies the 
first principle of thermodynamics. In the non-stationary case a volume 
^— «j . term appears since the solution does not admit a Killing vector. 

m ■ 

^ ' 1 Introduction 

O 

It is known that in many situations the entropy of a black hole solution 
^ ' of Einstein's equations can be calculated in the framework of a classical 

field theory without resorting to a statistical approach. 
This may be considered a good feature at present since statistical ap- 
proaches are based on a Hamiltonian and/or quantum formulation and 
, both these aspects are not yet clear for General Relativity. The definition 

» I ' of black hole entropy for stationary black holes as a macroscopical quan- 

I tity which satisfies the first principle of thermodynamics has in fact been 

proposed in [1], it has been settled on a secure mathematical ground in 
[2] and supported by a lot of examples [3], [4]. In this paper we aim to 
discuss a possible extension to non-stationary cases. 

We use the framework of classical field theory in its Lagrangian formula- 
tion; we require the theory to be natural, which means that each diff'eo- 
morphism of spacetime is a symmetry for the Lagrangian. For our purpose 
it is useful to rely on the geometric language of fiber bundles in which the 
calculus of variations is most naturally defined. In this framework, be- 
cause of the first principle of thermodynamics, the product between the 
"temperature of the black hole horizon" and the "variation of entropy" is 
equal to an integral of a suitable (n — 2)-form at space infinity (in a space 
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time of dimension n > 3). The integrated form a{L, a, X) is obtained 
as the variation of a suitable conserved quantity (in the sense of Nother's 
theorem) and it is associated to a vector ^ — 81 + ^^84, on spacetime. In the 
case of stationary black holes, it is possible to transform the integral at 
space infinity into an integral on a trapping surface 11 for the singularities, 
since the form a turns out to be closed: 

TSxS= [ a{L,a,^,X) (1) 

If the horizon of the black hole is bifurcate and we choose 11 to coincide 

with the bifurcation surface, we obtain, as a particular case, the same 
formula obtained by R. Wald and V. Iyer in [1]. We stress however that 
the above assumptions on 11 are unnecessary, often difficult to deal with 
and sometimes impossible (see ref.s [2], [3], [4]). 

In this paper we try to extend the same definition of entropy to the case 
of non-stationary black holes. We here consider black holes with an oscil- 
lating horizon without a quadrupole momentum, so that they do not emit 
gravitational waves (see, for example, R. Wald and V. Iyer [1] or Prolov 
[11], [12]), for which the first principle of black holes thermodynamics has 
the form: 

SxM = T5xS + mxJ (2) 
because they can be considered as isolated systems. 

The problems which occur in the non-stationary case are related to the 
fact that the vector ^ is no longer a Killing vector for the solution a, a 
fact that in stationary case is extensively used to prove that the form a is 
closed (see ref.s [1], [2]). When we try to transform the integral at spatial 
infinity into an integral on a " finite surface" , a volume integral of the so- 
called symplectic form uj{a,£^,X) appears because the form a is no longer 
closed, i.e. Div a{L, a, ^, X) ^ 0. In this case the entropy is defined as 
the sum of two integrals: 

TSxSdyn= [ a{L,a,i,X)+ [ Lo{a,^,X) (3) 

where E is the volume of the region enclosed between 11 and space infin- 
ity. We will show that both integrals axe well defined and can be easily 
evaluated in the framework we use. 

An earlier proposal for dynamic black hole entropy was given by R. Wald 
and V. Iyer [1]. They just tried to adapt the definition given for the sta- 
tionary case but, in the end, their definition turned out not to be covariant 
due to unessential requirements made on the existence and structure of 
the horizons (see [2]). The new definition we propose satisfies all the con- 
ditions imposed by R. Wald and V. Iyer in [1] and in addition, as a conse- 
quence of the geometrical framework we use, our proposal is automatically 
covariant with respect to any fibered morphism, i.e. any redefinition of 
fields. This is a stringent requirement from a physical viewpoint, since 
any lack of covariance produces results which are either wrong or at least 
require a lot of efforts in order to get a correct physical interpretation. 
A further problem, which we are not able to overcome, is due to the lack 
of examples on which one can test any prescription for entropy of non- 
stationary black holes. Basically we don't know any non-stationary and 
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geometrically well defined exact solution, for which one has a reasonable 
physical interpretation. For this reason we shall not present any example 
and direct application of our framework. Nevertheless we believe that the 
result we obtain is of interest, since it enlights the concept of entropy even 
in the stationary case by clarifHng which are the fundamental properties 
of entropy and which are instead mere consequences of stationarity. 
Furthermore, even if we do not know any explicit solution to test the for- 
malism, we stress that the class of solutions that are under consideration 
is certainly not empty, physically relevant and, as mentioned above, it has 
been taken into account in the literature (see references [1], [11], [16]). 

2 Notation and review of the stationary case 

Hereafter we recall briefly the standard notation and the definition of 
entropy in the the stationary case (more details can be found in [2], [5], 
[8], [9]). Let us consider a configuration bundle B fibered on a spacetime 
M and let us denote by j''{B) its fc-ordcr jet bundle, i.e. the space where 
fields live together with their derivatives up to order k included. Fibered 
local coordinates on j''{B) are defined by (a;^, J/', t/'p, y'^j^...,^^)- Let 
us also denote by A"(T*M) the bundle of n-forms over M. A Lagrangian 
of order k defined on B is a morphism of fiber bundles: 

L: J*(B) ^ A"(T*M) (4) 

In the case of General Relativity in vacuum we can choose the second 
order Hilbert Lagrangian: 

L = ^RVads (5) 

The variation of the generic Lagrangian (4) can be expressed through the 
so called first variation formula. Wc consider a vertical vector X on B 
and the variation of L along the flow of X, evaluated on a section a of B: 

< SL o j^alj^X >=< E{L) o falX > +d[< F{L,y) o f-^lj^-'X >] 

where E{L) and F{L,'y) are well-defined global morphisms. The Euler- 
Lagrange morphism: 

E{L) : J'"'{B) ^ A"{M)(g)V* (B) (6) 

is unique and it defines the field equations E(L) o j^''a = 0. In this 
case V*{B) denotes the dual bundle of the vector bundle V(B) of vertical 
vectors on B. The Poincare-Cartan morphism F{L, 7) depends in general 
on the Lagrangian and on an arbitrary background connection 7 on M: 

F{L, 7) : J^''-\B) A"-' (M) ® V* ( J^-^B) (7) 

In particular, the Eulcr-Lagrange morphism for General Relativity gives 
the vacuum Einstein equation 

< E{L)\X >= e^^.dg'^'ds = {R^. - \-^g^.)6g'''' ds (8) 
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while the expression of the Poincare-Cartan morphism is in this case 

<]F{L)\j'X>=PfgVeSg^''ds, (9) 

where we have set 

Kp = - (y^) V9 [9"' 90,0 - 5Ul)] (10) 

We can notice the Poincaxe-Cartan morphism (9) does not depend on any 
background connection due to the low order {k = 2) of the theory. 

We recall also that one can introduce a definition for the Lie derivative 
of bundle sections with respect to the flow of a vector H on B projectable 
over ^ on M as: 

i^HCr = Tct(0 - H o cr (11) 

A projectable vector field H is an infinitesimal symmetry of the Lagrangian 
iff the following holds: 

<5Lo I i''£s<T >= d{HL) (12) 

A bundle is natural iff for each spacetime diffeomorphism f : M M 
on the basis M it is possible to find a canonical lift 0/ : B — > B on the 
bundle. 

A field theory is natural iff the configuration bundle is natural and each 
diffeomorphism on the basis M is a symmetry of the Lagrangian L (in the 
sense that L is invariant under the pull back via any lift of 0/). In natural 
theories wc can define a Lie derivative with respect to a spacetime vector 
field ^ by setting: 

£^a = £b(t (13) 

In this letter we will treat only natural theories, but the formalism intro- 
duced here is also valid for the more general case of gauge-natural theories 
[2], [10]. 

In the case of the Hilbert Lagrangian (5) which gives to General Relar 
tivity the structure of a natural theory, the covariance condition (12) with 
respect to the vector ^ can be expressed in the form: 

d,{eL) = le^a^g'^" - ^V99"^£iRc.0 (14) 

where e,j^ are the coefficients of the Euler-Lagrange morphism given by 
equation (8). From the first- variation formula and the covariance condi- 
tion it is possible to formulate the Nother's theorem which associates to 
any vector field ^ on the spacetime, a conserved current £{L,^) so that: 

Div{£{L,^))=W{L,^) (15) 

where Div denotes the (formal) divergence operator. 

The work-form W{L, ^) vanishes on-shell, i.e. along solutions of field 

equations. For General Relativity in vacuum we obtain that: 



f'(i,0 



+ {Rp^-lV99p.)9''\'' (16) 
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For each natural theory (as well as for gauge-natural theories) using 
Bianchi's identities, it is possible to decompose the current £ as: 

E(L,i,)^l{L,(,)^Div{U(L,(,)) (17) 

where £ is defined the reduced current and U is defined the superpotential 
of the theory. The reduced current £ vanishes on-shell. For the Lagrangian 
(5) the gravitational superpotential U can be explicitely calculated (see 
[5]) and it is known to be: 

W(i,0 = -^Vs(V*^r-V,'^^^)dsv (18) 

which is called the Komar superpotential [5]. The conserved quantities 
can now be obtained integrating the current £ on a {n— l)-region D, i.e. 
a compact submanifold of M with a compact boundary dD. So the con- 
served quantities are integrals of the superpotential on dD. Generally the 
quantities obtained are not conserved with respect to the "time" defined 
by an ADM splitting of spacetirnc in spacelike surfaces. This happens, 
e.g., when the timelike vector ^ is a Killing vector for the solution g (see 
[18]). We stress that all the above quantities are linear with respect to 
the vector ^ (toghcther with its derivatives). 

However, if we calculate the conserved quantities for General Relativity 
integrating the superpotential (18) on dD, the mass obtained does not 
assume the physically expected value. This is the well known anomalous 
factor problem which affects the Komar superpotential. There are at least 

two different ways to solve the problem [5]. 

If we consider the variation of conserved quantities, defined above, this ex- 
pression will suggest us to define the variation of the corrected conserved 
quantities by the ADM prescription [5]: 

SxQD{L,^,a)^ [ [5xU{L,t<r)-id<F{L,l)°f''-'<r\j''-'X>)] 

J dD 

which gives us the expected quantities, of course up to an integration 
costaiit. To coiistmct this formula in a covariant way it is necessary to 
introduce a background connection 7 and the conserved quantities will de- 
pend on the background connection chosen. We can consider 7 to provide 
us a "zero level" for the energy, so that (as it is physically resonable) it 
is like a parameter for the theory. In the case of General Relativity we 
can choose as a sort of "natural" background connection the Levi-Civita 
connection of any background metric. 

On the other hand, one can notice that the Lagrangian (5) of General 

Relativity can be written as the sum of a first order Lagrangian and 
a divergence depending on the background metric. The first order La- 
grangian gives us a mechanism analogous to the ADM formalism (thought 
explicitely covariant and independent on the choice of a foliation) to cal- 
culate the corrected conserved quantities [6] . The background fixing pro- 
duces in this case an additional boundary term which solves the anomalus 
factor problem. The first order Lagrangian method is more general (in 
fact it is also applicable to non-compact solutions), while ADM formal- 
ism is not applicable in this case because space infinity is not a priori 
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asymptotically flat. This is particularly important to our purposes, since 
cosmological solutions are usually non-compact and have general a differ- 
ent asymptotical structure. 



Mass and angular momentum for a black hole solution of Einstein 
equations arc defined as the conserved quantities respectively connected 
to the vectors dt and 9^ on spacetime. The corrected quantities, which 
provide us the physically expected values, are defined as intcgals at space 
infinity: 

M = I [UKomar{L,dt,g) - B{L,dt,g) (19) 

J = - f [UKomar{L,d^,g)-B{L,d^,g)] (20) 

where the {n — 2)-form B is defined through integration of the variational 
equation on dD: 

5xB{L, t g) = i(Div < F{L, 7) o fg\j^X > (21) 

see for example [2], [7]. 

A covariant (and somehow canonical) choice for vacuum General Relativ- 
ity is: 

B{L, {, g) = -^g'-^i^'^w^dsx^ (22) 

where we have set 7 and V to be the ChistofFel's symbols for the metric 
and the background connection respectively and we have defined: 

f K, = <p - K, 

I f/M _ pM _ pP 

K ^al3 — Ial3 'p(f}°a) 



3 Definition of entropy in the stationary case 

The entropy of a stationary black hole solution is defined as the macro- 
scopical quantity which satisfies the first principle of thermodynamics (1). 
We impose that X is a solution of linearized field equations and both the 
temperature T and the angular velocity of black hole horizon Q are con- 
stant parameters depending on the class of solutions chosen. The tem- 
perature T is just defined as the temperature of the Hawking radiation 
T = where n is the surface gravity, as shown in [14], [18], [19], [20] 
by means of Euclidean path integrals. On the other hand Q is defined so 
that I ^ j'^ vanishes on the BH horizon. If we solve (1) with respect to 
SxS we obtain that: 

6xS=^{6M-mj) ==^ I l5xU{L,i,a)-ii<F{L,^)\fX>] 



where 00 means the space infinity of a spacelike slice and ^ = dt + ^d^, . 
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Under the hypotheses that ^ is a Kilhng vector for g and X is a solu- 
tion of hnearized field equations it is possible to prove in a very general 
framework (see [2]) that the quantity under integral: 

a{L,g,i,X) = 5xU{L,i,g) - < ]F{L,'~i)\j^X > (23) 

is a closed form. This allows us to redefine 5xS as an integral on any spa- 
tial surface which is homologically equivalent to cx). In this definition we 
do not have any additional requirement about maximality of the solution 
considered neither about the horizon properties. In particular it is not 
necessary to require ^ to vanish on the trapping surface. This fact allows 
us to apply the definition to a wider range of solutions and simplifies both 
conceptually and computationally the calculations (see [2], [3], [4]). 
If the solution admits a bifurcate Killing horizon and we can choose a 
bifurcation surface on which ^ vanishes, then our more general definition 
reproduces, as a very particular case, the one given by Wald and Iyer 
in [1]. This latter definition is not applicable to solutions for which ^ is 
not a Killing vector (non- stationary solutions, non-asymptotically flat 
solutions, etc.). 

4 Vairiation of conserved quantities 

It is possible to express a bundle morphism (for example the Euler- 

Lagrange and the Poincare-Cartan morphisms which are fe-forms on M) 
in local fibred coordinates. In this local formalism, any such morphism 
appears to be a linear combination of the vector field components (^'', 5') 
together with their derivatives up to order r (r = 1 for the example of 
General Relativity under investigation). Let us thence consider a deriva- 
tion 5 (i.e. a linear operator which satisfies the Leibniz rule). If we are 
able to calculate the J-derivative of a vector field component and of the 
whole n-form then, by applying the Leibniz rule, we are able to define the 
derivative of the coefficients of the linear combination. 
In our case we apply this rule to the Lie derivative and to the variation 
along the fiow of a vector field X, which are both derivations. For exam- 
ple, if we choose as a particular derivation onto the Poincare-Cartan 
morphism, where C is a vector field on spacetime and X € V{M), we 
obtain for a second order theory: 

£i < F{L,^)\fX >= £i[p^X' +p^''X;'\dsx (24) 

and applying the Leibniz rule we obtain: 

£iP7 = {d^ep7 - d^Cpt" + +ed^.p7 + di^pY - d^^pr) 

We remark that this definition for the Lie derivatives of the Poincare- 
Caxtan morphism can also be obtained by considering pt ,p7 as the local 
expressions of a section on a suitable fiber bundle and thence applying 
the general definition (11) of Lie derivatives of sections of fiber bundles 
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(see ref. [21], [22]). 



It is now possible to analyze the divergence of the form a{L, a, ^, X) 
in the case that no conditions whatsoever are imposed on the solution 
a. In the case of stationary black holes one requires ^ to be a Killing 
vector of the solution. Accordingly we assume in the general case that 
^ is a symmetry for a, i.e. £^a = which is fundamental to prove that 
da{L,(T,$,X) =0 (see [2]). 

If we relax this condition such a divergence does not vanish anylonger. 
This is the case of non-stationary black holes in a relativistic theory. In 
this general setting we have that: 

Div{a{L,^,a,X)) = Divi5xU{L,^) - < F{L,'y)\j''-''X >) = 

= 5x< F{L,y)\j''-'£^a > -£«(< F{L,j) | /"^X >) 
- Sx£{L,^)-i^<E{L)\X > 

(see for example [2]). 

Lot us thence analyze each term of this expression. The Euler-Lagrange 
morpliism vanishes on-shell < E{L)\X >= 0. The variation of the re- 
duced current can be expressed as: 

Div5x£{L, = -Sx < E{L)\£^a > (25) 

and this term is identically vanishing since X is a solution of the linearized 
field equation (see [2]). 

Using the prescription given for the variation of fiber bundle morphisms, 
in the case of a theory of order fc = 2, we can analyze the two terms left 
on the right hand side and we obtain as a special case: 

Div [a{L,^,a,X)]=cj{L,^,a,X) = 
= < 5xF{L)\j'^£ia > - < £iF{L,j)\j^X > + < F{L)\j^Z > 

where w is an (n — l)-form on spacetime and Z = Z^di is a vertical field 
defined as = {djX^£^a'^ so that it lifts to Z = Z^di + {d^,Z''^^^^^■ . 
Let us stress that in the above expression each term is "under control" in 
the sense that it can be analytically calculated whenever a Lagrangian is 

given for the theory. The expression of Div[a[L, ^, a, X)] is fundamental 
to our purpose; in fact it will contribute to the entropy formula under the 
form of a volume integral. 

To summarize, in the case of under analysis (General Relativity in vac- 
uum), we see that it is possible to calculate w and a, using the formula 
for the Komar supcrpotcntial and for the Poincarc-Cartan morphism, ex- 
pressed in local coordinates by (9), (18), namely in our case: 



f w(L,e,a,X) =5x< F{L,j)\j'£^g > -£^ < F{L,-i)\j^X > 
\ a{L, ^, a, X) = SxU{L, ^, g) - < F{L, 7) o fg\fX > 

where variations and Lie derivatives can also be defined according to (24) 
as usual for differential forms. 
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5 Non-stationciry black holes 



In this section wc extend our definition for black liolo entropy to the case 
of non-stationary black holes. As a motivation let us mention that cosmo- 
logical solutions of Einstein's equations are usually not stationary and not 
asymptotically flat models, i.e. the solution does not admit any timelike 
Killing vector. In this case the (n — 2)-form a is no longer closed and 
we cannot easily define the entropy as a boundary integral on a trapping 
surface for the singularity. 

In our model the first principle of thermodynamics is the same used for 
stationary black holes. We consider black holes which do not omit gravita- 
tional waves, so we consider only solutions without a quadrupole momen- 
tum [1], [11], [12]. This means that the system is isolated and electrically 
not charged. The geometrical formalism we use is manifestly covariant. 
We will also show that our proposal satisfies all the reliability conditions 
stated by Wald and Iyer in [1]. 

We define again the entropy for a non- stationary solution of Einstein equa- 
tions as the macroscopical quantity which satisfies the first principle of 
thermodynamics (1). The definition is the same given before for the case 
of stationary black holes. In this new case, however, T and Q cannot be 
calculated as the temperature and the angular velocity of black hole hori- 
zon, but they can be considered as a priori parameters of the theory. The 
only requirement is to ask these parameters to realize am integrable first 
principle of thermodynamics. The choice among them has to be carried 
over on the basis of some external physical consideration. However this 
is not a feature of non-stationary solutions; even in the stationary case, 
if we choose quasi-local energy instead of mass we obtain a different (but 
integrable) first principle [18]. This fact will be subject of further investi- 
gations. 

On the other hand, the mass and the angular momentum may no longer 
be time-conserved on a spacelike ADM slice of spacctime, but they are 
covariantly conserved, i.e. they obey a continuity equation; in other words 
they are conserved in the sense of Nother theorem even if their values may 

change in time. 

We can substitute in the first principle the expressions (19) and (20) for 
mass and angular momentum calculated by means of Nother theorem and 
we will obtain an expression which defines the variation of entropy as an 
integral on space infinity (23). Now, in the case of non-stationary black 
holes, we have to take into account that a is not closed to evaluate the 
same quantity on a trapping surface. We have to consider the form ui and 
its integral on a volume E between the trapping surface 11 and the space 
infinity. So we will obtain the formula for entropy under the form: 

T5xSdyn= [ a{L,g,0+ [ oj{L,g,^) (26) 

where T is the black hole temperature and the integrated forms w and a 
are defined in the previous section. In this formula each term is explicitly 
known and computable once the Lagrangian and the exact solution g are 
specified. The only restriction on the theories we analyze is the fact that 
they have to be well-defined from a Lagrangian viewpoint. 
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The effort to apply this formula to known solutions has boon vain up 
to now. In the case of General Relativity in vaeuum, in fact, to our 
knowledge there are no well defined exact non-stationary solutions in lit- 
erature, even if the Hilbert Lagrangian which defines the theory is global 
and covariant. If we consider otherwise the case of General Relativity 
in interaction with matter it is then possible to find in literature some 
explicit non-stationary solution; in this case it would be easily possible 
to generalize the definition of entropy (26) to treat also these theories by 
just adding am interaction term to the superpotential which eventually 
enters the final formula (26) (see [13]). In this latter case, however, there 
is no well defined global Lagrangian for the theory, because of the exhotic 
properties of the gas matter considered, which is an essential requirement 
to calculate conserved quantities in a geometrical framework. However 
we carry over a theoretical analysis of the case of General Realtivity in 
vacuum. 

Wald and Iyer imposed some conditions a priori on the reliability of 
the definition which our definition satisfies by default: 

• In the case of stationary black holes we must have T5xS = TSxSdyn- 
To show this fact it is enough to say that when the solution is stationary 
we have £^g = => ijj{L, g, ^) = 0. 

• We have to show that in the case of non-stationary perturbations, gen- 
erated by a field X, of a stationary solution we have TS^Sdyn = TS^S. 
In this case once again u}{L,g,£^) = 0. This is related to the fact that 
what we need is £(go = where go is the unperturbed solution, which is 
stationary. 

• The entropy for a theory defined by an equivalent Lagrangian L + DivO 

should be the same calculated for the theory defined by the Lagrangian 
L. It easy to see that for any pure divergence Lagrangian we have 
a{Div9,a,^) = 0. Since a is linear in L, in fact we have that a = => 
a; = because we have chosen a pure divergence Lagrangian. 

• Our definition must satisfy the second principle. This point is left out 
to future investigations, but we stress that at the moment the problem is 
still out of control even in the stationary ease [1] ; however it is reasonable 
to say that the second principle is related to the second variation of the 
fibered morphisms we have constructed, and thence to the positivity of 
energy. 

• Finally, our definition should be covariant under field redefinition. Our 
formalism is manifestly covariant by construction and our definition sat- 
isfies this condition too. It is easy to prove this claim if we consider the 
trasformation rules for the Poincare-Cartan morphism coefficients and for 
the field /"^X. 



6 Conclusions and perspectives 

We have proposed a new prescription to calculate the entropy for a non- 
stationary black hole. This formula is applicable to a well defined relativis- 

tic theory with a known (global) Lagrangian and whenever Lie derivatives 
are well-defined so that we can define covariant conserved quantities at 



10 



space infinity. The application of our formalism to calculate entropy for a 
well-defined explicit solution is immediate; once we have the solution in a 
local coordinate system it is possible to apply the algorithmical formalism 
we have developed to calculate the conserved quantities and entropy. The 
formula we have proposed is independent on the choice of the trapping 
horizon for the singularity; if wc consider horizons belonging to the same 
homotopy class, the result obtained for entropy is invariant (see for ex- 
ample the TAUB-BOLT solution, a discussion of which is given in [4]). 
Prom a physical viewpoint this definition satisfies the conditions stated 
by Wald and Iyer in [1] and in particular it is covariant. A future task 
will be to calculate explicitly the entropy for some exact non-stationary 
solution of Einstein equations. We will furthermore investigate the second 
principle for our definition both for stationary and non-stationary black 
holes. 
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